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I. INTRODUCTION 



Topological defects are important in condensed matter physics since they can affect prop- 
erties of matter or even the phase structure of the system. This can also be the case in the 
condensed matter physics of QCD jl|. It seems likely from theoretical studies that dense 
and cold quark matter is a color superconductor [2]. Such a state of matter can be realized 
in the core of neutron stars or possibly in heavy ion collisions. It is quite important to in- 
vestigate the properties of a color superconductor to verify its existence experimentally and 
observationally. In particular, at extremely high density and low temperature, perturbative 
calculations are reliable, and it is believed that color-flavor locked (CFL) phase emerges 
where all of the three light quarks make pairs and form condensates 0, y|. The original 
symmetry of QCD, G = U(1)b x SU(3) c x SU(3)l x SU(3)r, is spontaneously broken down 
to the color-flavor locked symmetry H = SU(3) c+ l+r = SU(3) c+ p in the CFL phase, apart 
from discrete symmetry. It exhibits superfluidity due to the breaking of the global U(1)b 
symmetry, as well as color superconductivity because of broken color symmetry. 

By examining the topology of order parameter space G/H ~ U(3) it was shown for the 
first time in Ref. ff| that topologically stable vortices exist as a consequence of the symmetry 
breaking in the CFL phase. The topologically stable vortices with the lowest winding number 



in the CFL phase are non-Abelian semisuperfluid vortices. They are ca 
in the sense that the remaining symmetry if is a non-Abelian group 



led "non-Abelian" 
1. Just vortex 



in a conventional type II superconductor carries a magnetic flux, a non-Abelian vortex in 
the CFL phase carries a color magnetic flux. At the same time, this object behaves as a 
superfluid vortex as a result of superfluid properties of the color superconductor. Therefore, 
when a color superconductor rotates faster than some critical velocity, these vortices are 
created along the axis of the rotation. When the rotation speed is further increased, the 
created multiple vortices should form a lattice of vortices as in ultra cold atomic superfluids. 
As a result, if the density in the core of a neutron star is sufficiently high for the color 
superconducting matter to exist, a vortex lattice should appear and it should yield some 
physical consequences which can give a signal. The superfluid turbulence may also occur if 
there is some mechanism which leads to an instability of ordered structure like the Donnelly- 
Glaberson instability j§| in atomic superfluids. 

In order to find out what kind of collective structure of vortices is realized, it is essential 
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to determine the interaction between vortices. This can be clarified by investigating the 
interaction of vortices with quasiparticles, that mediate the force between vortices. The 
interaction between vortices and quasiparticles is also useful to study the transport properties 
of a color superconductor. As quasiparticles, there appear eight gluons with masses given 
by the Higgs mechanism, eight Nambu-Goldstone (NG) bosons (CFL mesons) associated 
with the chiral symmetry breaking, and one NG boson (phonon) due to the breaking of the 
baryon number symmetry U(1)b- 

One striking feature of non-Abelian semisuperfluid vortices is that they have internal 
degrees of freedom, which are called orientational zero modes. The existence of a vortex 
breaks the color-flavor locked symmetry SU(3) c+ f down to its subgroup SU(2) x U(l) around 
the core of the vortex. Consequently, there appear further NG modes confined inside the 
core of the non-Abelian vortex [9|, which parametrize the coset space 

2 _ SU(3) C+F (1) 

" SU(2) xU(l)' 1 } 

There are degenerate vortex solutions with different color magnetic fluxes, which correspond 
to points on the CP 2 space. The force between largely separated non-Abelian vortices in the 
CFL phase is independent of internal orientations [9] . This can be understood by noting that 
the long-range force is mediated only by massless U(1)b phonons, which are insensitive to 
color fluxes. However when the separation between vortices is relatively small, one expects 
that the force can be mediated also by massive particles which are not color singlet such 
as gluons, and consequently the short-range force can depend on the color fluxes or the 
orientational zero modes of the vortices. 

In this paper, we determine the interaction between a non-Abelian vortex and quasipar- 
ticles such as gluons and phonons. To this end we derive a dual Lagrangian corresponding 
to the low-energy effective theory of the CFL phase by making use of a dual transformation. 
A dual transformation is a method widely used in various fields of physics, which relates 
theories that have different Lagrangians and variables with the roles of equations of motion 
and Bianchi identities exchanged. Under this transformation, particles and solitons typically 
interchange their roles, i.e. Noether currents are interchanged with topological currents, and 
vice- versa. Topological defects in the original theory are mapped to fundamental particles in 
the dually transformed description. One of the advantages of dualization is that one can deal 
with the interaction of topological defects by the method of ordinary field theory for parti- 
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cles. Dualization was first used to deal with Abelian vortices; in the case of local vortices in 



the Abelian Higgs model, the U(l) gauge field is mapped to massive two-form (antisymmet- 



ric tensor) fields 



10|. On the other hand, in the case o 



U(l) NG mode is mapped to massless two-form fields 
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obal (or superfluid) vortices, the 
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mapped to strings which behave as a source of two-form fields 



n either case, vortices are 



ll| . Duality can be extended 



to non- Abelian gauge field theories; The dual transformation was applied by Seo, Okawa 



and Sugamoto 



13| to SU(2) gauge theory of Higgs fields in the adjoint representation, which 



admits the so-called Z 2 vortices. Gauge fields are dualized to non-Abelian generalization 



of two-form fields 
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14j and the Z 2 vortex appears as a source of non-Abelian two-form 



fields. However the Z 2 vortices are Abelian and do not have internal orientational modes, 
although the gauge fields and dualized two-form fields are non-Abelian. 

We apply the duality to non-Abelian vortices in the CFL phase, in which internal orienta- 
tional modes reside. Whether or how these orientational zero modes couple to quasiparticles 
is a nontrivial question. We show that the interaction of non-Abelian vortices with quasi- 
particles naturally arises by a dual transformation. In the dual theory, dual gluon fields 
are described by massive non-Abelian two- form fields as in 13]. We find that non-Abelian 
two-form fields are coupled to non-Abelian vortices through their internal orientational zero 
modes given in Eq. ([I]). This result is quite natural because these modes correspond one-to- 
one to the color magnetic flux which the non-Abelian vortex carries. On the other hand, the 
phonons are dual to Abelian two-form fields coupled to non-Abelian vortices through their 
translational zero modes in the same way as for Abelian vortices, but the coupling strength 
is 1/3 of that for Abelian vortices. Because the interaction terms do not involve the space- 
time metric, they are called topological interactions. The dual Lagrangian obtained here 
provides a starting point to analyze a possible collective structure of vortices. 

This paper is organized as follows. In Sec. [Ill we give a Ginzburg-Landau Lagrangian 
of the CFL phase. In Sec. II I II we perform a dual transformation of the Ginzburg-Landau 
Lagrangian. First we take a dual of massive gluon fields to massive non-Abelian two-form 
fields. Then we take a dual of the U(1)b phonon to the massless two-form field. In Sec. IIVI 
we use a solution of a single non-Abelian vortex to calculate the vorticity tensor of the 
non-Abelian vortex explicitly. Sec. |V] is devoted to summary and discussion. In Appendix 
we will give some details of calculations omitted in the main text. 
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II. LOW-ENERGY EFFECTIVE THEORY OF THE CFL PHASE 



In this section we introduce an effective theory of the CFL phase and discuss the sym- 
metry of the ground state. We start with a time-dependent Ginzburg-Landau(GL) effective 
Lagrangian |29| for the CFL phase up to second order in time and spatial derivatives. The 
Lagrangian is written in terms of order parameters <3> L and $ R , 



L(R)l 



L(R)^-y L(R)\ 



eabctijk \q bj Cq ck J 



(2) 



where a, b, c and i, j, k are color and flavor indices, respectively, and g L ( R ) are quark fields of 
left (right) chirality. The crossing terms of <£> L and $ R are allowed by symmetries, but are 
suppressed at high densities Here we take $ L = — $ R = $ so that the ground state is 
positive parity state. In the following we neglect the effect of U(1)em electromagnetism since 
the mixing between broken SU(3) C color and U(1)em is sufficiently small at high densities. 
We also consider a sufficiently high density region where the masses of three light quarks can 
be neglected. (We refer the reader to |20] for a recent discussion on the effect of a strange 
quark mass.) 



In this case the GL Lagrangian for the CFL phase is given by 
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17| 



C(x) = l(E a ) 2 -^(B a ) 2 + ^oTr [(D $) + D°$ 
V($) = Tr [A 1 ($ t $) 2 -A 2 <I> t <l>] +A 3 (Tr[$ t $]) 2 , 



Ki Tr 



(A$) t + iK' Q Tr [$ f £>o$] 



(3) 



where Ef = F& Bf = \e ijk F« k , F* v = d,A a u - d u A* + gf abc A b ^ , = (d, - igA«T a ) $ 
(i — 1, • • ■ ,3), and T a are generators of SU(iV) c normalized as Tr[T a T b ] = \5 ab with color 
indices a = 1, 2, • • • , N 2 — 1. Here the Lagrangian describes the low-energy effective theory 
in the CFL phase only for N = 3, but we have extended the order parameter field $ to an 
N by N matrix. Coefficients K , K' , Ki, X%, A 2 , and A3 are GL parameters dependent on 
the temperature and the chemical potential of the system but are dealt with as constant 
parameters in this paper, e and A are the dielectric constant and the magnetic permeability. 
The form of the kinetic term of gauge fields in the effective Lagrangian is deduced by 
requiring the gauge invariance, rotational invariance and parity conservation. The Lorentz 
symmetry does not have to be maintained in general since superconducting matter exists. 
However, there exists a modified Lorentz symmetry in which the speed of light is replaced 



by 1/veA. It is always possible to restore the Lorentz invariance of the kinetic term of gauge 
fields by rescaling x°, Aq, K , K' q and K\. Therefore we can start with the Lagrangian in 
which e and A are taken to be unity. 

For notational convenience, we introduce a vector = (K , Ki, Ki, Ki). Our starting 
point is the following GL Lagrangian 

C(x) = —(F^Y + K^Tr^D^D^] +iK' Tr[^D ^\-V(^). (4) 

We consider the parameter region Ai > 0, A 2 > 0, Ai + A^A 3 > 0. This Lagrangian includes 
the most general terms which are consistent with the symmetry group G of generalized QCD 
with N colors and N flavors, 

Q = SU(iV)c x SU(iV) F x U(1) B 

(^Jv)c+B x (^Jv)f+B 

where SU(iV) c is the local color symmetry, SU(A r )p is the global flavor symmetry, and U(1)b 
is the global symmetry associated with the baryon number conservation. Under the action 
of the element (V c , V F) e l6 ) e G, $ transforms as 

$ _>. $' = e ie V c <Z>Vj, (6) 

where V c G SU(A^) C , V F e SU(N) F ,e ie e U(1) B . The elements corresponding to the discrete 
groups in the denominator of G can be written as (z±, Z2, Z\ 1 z 2 ^ 1 ) £ G with 21,^2 £ ^n- 
These transformations are removed since they do not change $ for any value of <3>. 

In the CFL phase the free energy is minimized when $ is proportional to a constant 
unitary matrix. By using the symmetry G, one can take the value of $ without loss of 
generality as 

$ = |A|1 W , (7) 

where |A| is a real number. By this expectation value, $ is invariant under the restricted 
transformations 

{h\h = (Uz~\ U*, z),U e SU(iV), zeZ N }c G. (8) 
Therefore the symmetry G is spontaneously broken in the ground state to 

H = (9, 

The order parameter space which characterizes the degenerate ground states is given by 

G/flg U(D^U(AQ 



The CFL phase admits stable vortices since the first homotopy group of the order parameter 
space is nontrivial: 

7Ti ~ Z. (11) 



III. THE DUAL TRANSFORMATION 

In this section, we perform a dual transformation within path integral to derive a dual 
Lagrangian for the CFL phase. After the transformation, massive gluons are described by 



massive non-Abelian antisymmetric tensor fields 13] and U(1)b phonons are described by 
massless antisymmetric tensor fields. We show that in the dual description vortices appear 
as sources which can absorb or emit these particles. This is consistent with the empirical 
rule that a dual transformation interchanges the role of particles and solitons. 

A. The dual transformation of massive gluons 

The partition function of the CFL phase can be written as 

Z = J VAl{x)V$(x) exp |i J d 4 x£(x)| , (12) 

with the Lagrangian defined in Eq. (111) . We shall impose the gauge fixing condition on the 
field $ rather than on the gauge fields since they are integrated out in the end. The gauge 
fixing condition is taken care of when we consider a concrete vortex solution. 

Let us introduce non-Abelian antisymmetric tensor fields B® u by a Hubbard- Stratonovich 
transformation 



exp 



J d 4 x|-i(F;j 2 | cc I VB; u e W i J d 4 * j-i [m 2 (B; u ) 2 -2mB; u F a n]] 



3) 



where = ^e^ upa B a ' pa . The parameter m introduced above is a free parameter at this 
stage. We will choose m later so that the kinetic term of B® is canonically normalized. 

Substituting ( IT3l) into (I12j) . we can now perform the integration over the gauge fields A^. 
The degrees of freedom of gluons are expressed by B a after this transformation. Each term 
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in the Lagrangian is transformed as follows: 



Ti{(D^Y(D^)} + iK' Tr [&D <f>] 
= Tr [^(V^ + igA a ^T a )(^^ - igA b '»T b )<&} + iK' Tr [<^> + (^ - igA*T a )<S>] 
= Tr{(9 M $) t (^$)} + iK' Tr [$ f 5b$] + gA^J^ 

+ff 2 VV^ Aa ^ MTr [$tT a T 6 $] , (14) 



with 31 = -iKp Tr 



$ f ( - 7^)T $1 - ^ Tr [$tT a $] , and 



1 



1 



-mB; u F a ^ = --mB;„(2d„A« + gf^A^) 



1 



(15) 



Performing the integration over A a , the following part of the partition function is rewritten 

as 

j PA^expji y d 4 x ^g 2 A a ^K^A h ' v - m (d v B« v - £ J°) A a '" | 
oc (det^-^exp U f ^ \-\ {^j (d p B^ - (K^)* (d a B^ - M-J>* 



(16) 



where K ab u is defined by 



K* = ^y/IQ^Tr [&T a T b <5>] - -f abc B 
^ 9 



c 

\1V 



jfyCib _ p>ab 



(17) 



with $J = y^y/K^T* [$ f T a T fc $] and Bf v = ! ahc B%. We define the inverse of K ab v 
by the power-series expansion in 1/g 

-1 oc 



X- 1 



9 



n=0 



(18) 



As a result, we obtain the following partition function 

Zoc| PflJUdetffJr^exp ji / d4 ^G^)}> 
where C* G denotes the gluonic part of the dual Lagrangian 



(19) 



(20) 
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Now we define the non-Abelian vorticity tensor uj^ v as the coefficient of the term lin- 
early proportional to . Collecting relevant terms in the above Lagrangian, the coupling 
between massive gluons and the vorticity is given by 



1 m 
27 

1 / m 
~2 17 



r>a , a, Act 
#Act W 





( m\ 




Kg J 



r^[$>- l B$- l ]f u j b < u 



(21) 



where we have defined the vorticity tensor u: a as 

,a,ACT = £ Act^ ^ {($-l)(J)jV} + je,a^-l)« yroda^-lj* j6,/9] _ 



(22) 



Here A(^ v ) is a symmetrized summation defined by A( ab v ) = A ab v + A ba . This expression 
for the non-Abelian vorticity is valid for general vortex configurations. The information of 
vortex configuration is included in <fr and J* 



B. The dual transformation of U(1)b phonons 

In the following, we perform a dual transformation of the NG boson associated with the 
breaking of U(1)b symmetry. This mode corresponds to the fluctuation of the overall phase 
of $ which can be parametrized as <&(x) = e l7T( - x ' ) ip(x), where ir(x) is a real scalar field. 



Substituting this into the following part in the Lagrangian (j4j leads to [30] 



K lt r Fr{{d lt $)\&>$)} + iK' Tr{&d 3>} = {d^f M 2 - 9"ttJJ + i^Tr(^) 2 + iK' Tr{i(>%i(>}, 

(23) 

with J° = S^oK' Q M 2 and M 2 = Tr [ip^ip] ■ We will transform the U(1)b phonon field tt(x) 
into a massless two-form field BP. Note that the field ir(x) has a multivalued part in general; 
since n(x) is the phase degree of freedom, n(x) can be multivalued without violating the 
single- valuedness of <J>(x). In fact the multivalued part of tc(x) corresponds to a vortex. Let 
us denote the multivalued part of ir(x) as tt mv (x). 

The dual transformation of this U(1)b phonon field is essentially the same as the case of 
a superfluid. We basically follow the argument of Let us introduce an auxiliary field 
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C M by linearizing the kinetic term of tt(x) in the partition function as follows 



Z oc / V-nVnMV exp % 



J d 4 x (M 2 K, t {d^TT + ir MV )} 2 - 0"(tt + ttmv) 4) 



oc / VtxVkmv'DC^ exp i 



r 2 

d A x 1 



(24) 



Integration over 7r(x) gives a delta function 



Jvirexpi y"d 4 x(-2C M ^:^7r + 7r^Jj) = 5 {d" (2C mV ^ + Jj) } . (25) 



Then let us introduce the dual antisymmetric tensor field B° by 



VC,VB%5 (2C,^K; + Jj - m°d»B% 



(26) 

where the dots denote the rest of the integrand and m° is a parameter. By this change of 
variables we have introduced an infinite gauge volume, corresponding to the transformation 
5B® U = d^A u — dyh^ with a massless vector field A /t . This can be taken care of by fixing 
the gauge later. There is no nontrivial Jacobian factor as the change of variables is linear. 
Integrating over C^, and transforming a resultant term in the Lagrangian as 



m°d v Bl^7r MV 



-m 



B^e, vp(T d v d^MV 



= -2nm a B°> pa uj pa , 
where the first equality holds up to a total derivative and we have defined 



(27) 



1 



d»d^ MV . 



pu 2tt 

We thus obtain the dual Lagrangian for the U(1)b phonon part 



(28) 



pii 



2M 



KJm°d v Bl - jy - 2nm»B^u;» 



t0\2 



nO,p, ,0 



(29) 



Note that the term linear in BP coming from the first term of (|2"9"|) is a total derivative and 
does not contribute to the equation of motion. The partition function is proportional to 



Z oc / Vn MV VB° u expi 



d A xC 



(30) 



The U(1)b phonons are now described by a massless two-form field B Q and vortices appear 



as sources for B® u . 
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C. The dual Lagrangian 



We now summarize the results obtained so far. We have shown that the partition function 
Z of the CFL phase is proportional to Z* with the dual Lagrangian C*\ 

Z oc Z* = J VB a ^Vn MV VB° u Vij (det Kj)~ 1/2 exp jz J d 4 x£*(x) J , (31) 

where 

C* = C* G + £* h + K, Tr(9^) 2 + iK' Tr{^d ^} - Vty). (32) 

Here C* G and £ Ph are given in (|2"U|) and (|2T?j) . respectively. We can discuss the interaction 
between vortices and quasiparticles in terms of the dual Lagrangian. Vortices are expected 
to appear as a source term for gluons and U(1)b phonons. The result above is valid for 
general vortex configurations. 



IV. COUPLING OF NON-ABELIAN VORTEX WITH DUAL FIELDS 

In the last section we have obtained the dual Lagrangian for general vortex configurations. 
In order to see the coupling of quasiparticles with internal orientational degrees of freedom 
explicitly and to discuss physical consequences, let us consider a single-vortex solution and 
find an expression of the vorticity tensor for this case. 



A. Vortex solution 



We consider a cylindrically symmetric vortex configuration along the z-axis. Then, the 



profile of a vortex solution with the lowest energy takes the form 



W=K ° |, (33) 
\^ <7(r)l*_ ly 

where f(r) and g(r) are functions of the radial coordinate which can be determined by 
solving equations of motions under their asymptotic behaviors (f,g) = (|A|, |A|) as r — > oo 
and (f,g') = (0,0) as r ->■ 0. 

Let us identify the CP N_1 zero mode in the background solution. To this end, it is 
convenient to take the singular gauge. We perform a gauge transformation on this solution 
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by V G SU(JV) C , 

( e -i6( 

V = I 

e^l^x 



-ie(N-i)/N n 
V=l " h (34) 



which transforms $ as 

$* = U$ 



f(r) 

g{r)l N -!) W) 

N-l 



= e ie / N { F{r)l N + G(r^ ' v 



iljv-i 



where F(r) and G(r) are functions of the radial coordinate which are related to / (r) and 
g{r) as 

f = F-^±G, g = F+±G. (36) 

Under this gauge transformation, the vortex solution is physically unchanged. At this stage 
we fix the gauge by fixing all the local color transformation of $ |19| . 

General vortex solutions are obtained by acting color-flavor locked SU(iV) c+ F transfor- 
mations on $*. Because of the existence of the vortex, the symmetry further breaks down 
to SU(iV — 1) x U(l). We parametrize the associated CP^ -1 orientational moduli space, 
which is obtained by performing the SU(iV) c+ F transformation on $* as 

$ = U^U- 1 = B ie ' N \f{t)1 n + G(r) Up - il^ | , (37) 

where U G SU(iV), is a complex TV-component vector which transforms as the fundamental 
representation of SU(iV) c+F and satisfies the relation 0*0 = 1. 

The definition of has a redundancy in the overall phase. We cannot distinguish and 
e"*0 because both express the same solution, so they should identified: ~ e"*0. Therefore 
are indeed the homogeneous coordinates on CP^ -1 . The low-energy excitation on the 
non-Abelian vortex can be described by the CP^ -1 model. It has been shown that the 
(Qpjv-i moc i es are j n f ac t normalizable and localized around the vortex 



19 



20|. 



We shall express the vorticity tensor by orientational zero modes and the profile func- 
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tions. can be written as 



(38) 



= -fig^y/K^l (V 6 + 2 ^T a T b ^j , 
where we have defined C(r) and D(r) by 

C(r) = (F — G/N) 2 , D(r) = G 2 + 2G{F-G/N), (39) 
whose asymptotic behaviors are (C,D) ~ (|A| 2 ,0) as r — > oo and 

(C,£>) ~ {(1 - 2/iV) 2 |A| 2 , (1 - 2/iV) 2 |A| 2 } as r ->■ 0. We define the inverse of $J by an 
expansion in the dimensionless quantity 2D/C. The latter is a deviation of $ from the 
ground state value and small except for the center of the vortex. Although we expand 
($ -1 )^ in power-series, we will sum up all the terms later. Therefore ($ _1 )^t is calculated 
exactly whenever the condition 2D/C < 1 is satisfied. This condition holds except for the 
vicinity of the vortex core. In matrix notation, $ reads 

* = (l + 2 4A > (40) 



4 V C 

where the components of the matrix and K are given by <f) ab = ^T a T b (J) and K^ v = 
g^v^/K^Ky . Then {&~ l ) ab v can be written as 

n=0 



More explicitly, 



— ^ (V - ^Vt°TV + ^T a T c <P^T c T b <f) + •••). (42) 

v K^K V c y C \ C J J 

On the other hand, the current J" is written as 

j; = -2^ImTr [$^7°$] - 5 m0 ^o Tr{$ f T a $} 

= -Ai^^ + ^^D(r)0t T ^ (43) 

B. The vorticity tensor 



Let us see how the vorticity tensor can be expressed in terms of profile functions and 
orientational zero modes in the case of a single-vortex configuration discussed in the last 
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section. The Abelian component ur of the vorticity tensor is readily identified as 

1 

2nN 



< = ^W^^W. (44) 



The non-Abelian component u;* of the vorticity tensor is obtained by substituting J® 
and into f[2"2"j) . In order to simplify the following calculations, let us define [l9j 

JF^a, b) = a<pd^ + bd^0 + (a - 6)# f 9 M # + , a, 6 G C. (45) 

This quantity J 7 satisfies the following relations 

b)^ = 7,{b\ a*), Tr[.F>, 6)] = 0, (46) 
aJ^(a, b) = J" M (aa,afe), ^(a, 6) + J>(a', 6') = J* M (a + a', b + 6'), (47) 
0+7>(a, 6)0 = 0, 0t [jr M ( a , 6), T a ] = Tr (^(a, -6)T a ) , (48) 
Tr{T a [^(a,6)t,J- H ( a ,6)]} (49) 
= (|a| 2 + |6| 2 ) {0t T a 0^0t^0 + d [u ^T a <p^ + ^T a d [u <P^d, ]( p + ^d [u ct>^T a d^} , 

(50) 

where [...] denotes antisymmetrization of indices. The current for a — 1, ... , N 2 — 1 can 
be written in terms of T as 

j; = —DK^ \d,9 + <^T> - 2G 2 K, ImTr {^T a d^)} 

= —DK^ {d,6 + ^|^o} ^T a <P + iG 2 K, Tr {T»J>(-1, 1)} ( 51 ) 

= W + W 

As is shown in the appendix, the following quantities can be written in terms of J 7 as 



"'■OS = -|70 f T a {^0 + g^ } + <a Tr [T^(-l, 



($- 1 J 1 )« + ($- 1 J 2 )«, (52) 



(J*- 1 )^ = -^ 7 f T a {^0 + ^<V>} - Mr [T^(/3, -1)] , (53) 

where the functions a(r), /3(r) and 7(r) are defined by 

AG 2 C D 

a(r) = —, m=c + n. M S C + D{1-W (54) 
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Hence the following part in the first term of the vorticity tensor reads 

+ = -^ 7 0t T V |a M e + ^%o} + «'a(l + 0) Tr [T°J>(-1, 1)] , (55) 

where we have used the linearity of J 7 . The second term in (|22|) is also expressed by J 7 , 
which is proportional to 



(Jcfr- 1 )^^- 1 J)5 = ((Jj + Ja)* -1 )^*- 1 ^! + Ja))5 
Each term is calculated as follows: 



(56) 



;ji*-V*-v 2 )^ + (Ja*- 1 /*" 1 ^)^ 

= -| 7 ^TV I V + ^^o} / abc Tr(T c ^(-a, a/3)) + ■ ■ 
= -07(1 + 0) ld,9 + -^5, j Tr {^(1, l)T a } , 



(57) 



(j 3 #- 1 )!;/ afc (*- 1 j a )» 

= - Tr [r 6 ^(-a/3,a)] / afcc Tr [T c ^(-«, «/3)] 



1 CM 

-— Tr{T a [^(-l,/3)t,^(-l,/3)]} 



(58) 



where we have used f abc T c = — i[T a ,T b ] and traceless property of J 7 . This quantity is 
explicitly rewritten in terms of as 

= (1 + (3 2 ) l^T a <pd { ^d v] <P + d^T a d v] 4> + ^T a d { ^d v ^4> + %0 f T a 0^]0V] , (59) 

where we have antisymmetrized the quantity with respect to (/i, v), because symmetric part 
vanishes when contracted with a completely antisymmetric tensor e\ atiU . Therefore non- 
Abelian components of the vorticity tensor is written in terms of J 7 as 



d v {-^Vt> (d»6 + + H 1 + 0) Tr [T^i-l, 1)] 

+ 1« 7 (1 + /?) ^ + Tr [TV^l, 1)] 



_!^Tr{T Q [^(-l,/3)t,^(-l,/3)]} 



(60) 
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This can be represented explicitly by orientational zero modes as 



"la = ^ 



+ia(l + (3) (fiT'd^cf) - d ti $T a (j) + 20 t T a 
-■^07(1 + 0) (d^^T a 4> + 4>^T a d^4>) (d u] 6 + ^^$ u]0 

-^a 2 (l + (3 2 ) (^T a (f)d^^d u] 4> + d [tl ^T a d v] <P + $T a d [ll <j)d v] $(j) + d [ »^T a 0d u] ^4>) 

(61) 

Now we estimate the relative importance of terms in (IBTj) . Let us define a parameter 
e = D/C. Since e expresses the deviation of the profile function from the ground state 
value, e is small away from the vortex core, as noted above. We can estimate the strength 
of each term of the vorticity in terms of the order of e: a(r) ~ 0(e 2 ), f3(r) ~ O(e ) and 
7(r) ~ 0(e). Therefore the leading order part of the vorticity is given by 



Acr 



+ 0(e 2 ). 



(62) 



In summary, the interaction of vortices and quasiparticles, dual phonons and dual gluons, 
is described by the action 

Sint — ^int + ^hrt, (63) 

respectively, in which <S?^ and S® t are defined by 

-2nm° I iV,r D\]„J l " !/ . (64) 



°int 



int 



m 
9 



d 4 x B a ^u a '^. 



(65) 



Here uj q and io a are given by (jHJ) and (IBTj) . respectively. Since these interaction terms do 
not couple to the space-time metric, they are topological interactions. 

The parameters m and m° are chosen so that the kinetic term of two-form fields are 
canonically normalized. We require the kinetic terms of B® u and BP to take the form 



d 4 x 



12K 



fiua 



(66) 



where H% a = d^ + d^ + d^, H% a = d^ + d^ + d^ and K, va = e p , V(T K" 

Under this requirement, the parameters m and m° are determined to be 

gVC{r)_ m o = ^__ ,/NC(r) + D(r) 
~ V2 ~ 



m 



(67) 
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The interaction is localized around the vortex while gluons and phonons propagate in 
the bulk. This means that the vortex actually appears as a source for these particles. Each 
term in (IBTj) is proportional to a(r) or j(r) which are functions of the radial coordinate and 
is nonzero only in the vicinity of the vortex. The functions a(r) or j(r) decay exponentially 
away from the center of the vortex with the characteristic length, minfmq 1 , m" 1 ), where htlq 



is the mass of gluons and m x is the mass of traceless part of the scalar field $ [18j . Specific 
form of the profile functions at(r), fl(r) and j(r) can be determined by solving the equations 
of motions numerically. 

Now we discuss some properties of the interaction (163 p . First let us look at the interaction 
(l64p of vortices with U(1)b phonons. This part is essentially the same as a vortex in a 
superfluid. For general vortex configurations we can write the Abelian component of the 
vorticity tensor as 

(a; )"" (s) = i J drd / ( ^'*V >(s - X"(t, a)), (68) 

where the the space-time position of the vortex world sheet is parametrized as X^(t, a) with 
the world-sheet coordinates r, a. The interaction of vortices with U(1)b phonons is written 

as 

^ = f d^B° (69) 



N 



where dcx Atiy = a , ' t - drdtr is an area element of the vortex world sheet. The magnitude 
of the Abelian vorticity is proportional to the winding number with respect to U(1)b, so 
we have the factor 1/N which is the winding number of vortices of the lowest energy in 
the CFL phase. U(1)b phonons BP do not couple to the orientational zero modes on the 
vortices, which is anticipated since BP is a singlet under SU(iV) c+ F while the orientational 
zero modes are in the fundamental representation. 

Next, let us look at the interaction of vortices with gluons ( )65|) . As can be seen in (|6"T1) . 
gluons actually couple to the orientational zero modes on the vortex. As a result, virtual 
gluons can be emitted and absorbed from orientational zero modes |3J|. This will lead to 
physical effects. For example, when two vortices are in a small distance, the orientational zero 
modes residing on them can exchange virtual gluons, which results in orientation-dependent 
corrections to the vortex-vortex interaction. It is also possible that the orientational zero 
modes on a single vortex emit and absorb virtual gluons. In the case of Abelian vortices, it 
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is known that the emission and ab sorp tion of virtual photons from a single vortex result in 
the renormalization of its tension [2l|. The emission and absorption of virtual gluons can 
lead to similar effects for non-Abelian vortices. 

V. SUMMARY AND DISCUSSION 

We have derived a dual Lagrangian of the Ginzburg-Landau effective Lagrangian for the 
CFL phase, in which massive gluon fields have been dualized to massive non-Abelian two- 
form fields while the U(1)b phonon field has been dualized to the massless Abelian two-form 
field. In the dual theory, non-Abelian vortices have appeared as sources of these dual fields, 
and we have thus obtained the topological interactions of non-Abelian vortices with these 
quasiparticles. By making use of a single non-Abelian vortex solution, we have explicitly 
calculated the vorticity tensor of the single non-Abelian vortex, and then have found that 
the phonons couple to the translational zero modes of the vortices as in the Abelian vortices 
while the gluons couple to the internal orientational zero modes CP N ~ 1 of them. 

Let us make comments on some applications of the dual Lagrangian obtained above. 

We can investigate the nature of the force between vortices. When vortices are separated 
at a large distance, U(1)b phonons mediate the most part of the force between them; the force 
between two vortices at the distance R is proportional to 1/R which can be understood 
as follows. The tension of individual vortex with the U(1)b winding k is E ~ A; 2 log A with 
the system size A. Therefore the tension of a non-Abelian vortex is 1/A 2 of that of a U(1)b 
vortex. Consequently the interaction energy for two non-Abelian vortices and likewise the 
force between them is 1/A 2 of those between two Abelian U(1)b vortices. This is consistent 
with our result f )69|) of the coupling of U(1)b phonons to non-Abelian vortices, which is 1/N 
of the one of U(1)b phonons to Abelian vortices; the force between two non-Abelian vortices 
is proportional to the product of the U(1)b charges of them, which is 1/N 2 . 

On the other hand, when the separation of vortices is small, exchange of virtual gluons 
will also contribue to the force between vortices. In that case the property of the force 
is dependent on the CP^ -1 internal orientations <p of the vortices because gluons (non- 
Abelian two-form fields) couple to the internal orientations through (IBTj) . This is natural 
since non-Abelian vortices carry color magnetic fluxes. 

It is also interesting to apply the dual Lagrangian we have obtained to the physics of 
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the core of dense stars. It is well known that the rotation of ultra cold atomic superfluids 
leads to a vortex lattice structure. Neutron stars are also rotating rapidly, so we can expect 
that a lattice structure of non-Abelian vortices will be formed if the core is dense enough 
to accommodate the CFL matter. If the speed of rotation is high, the distance between 
neighboring vortices is small and internal orientation will affect the nature of the interaction 
of vortices. As the type of collective structures of vortices is sensitive to details of the 
interaction between them, it is interesting to investigate the phase diagram as a function 
of the speed of rotations by making use of the action we have derived. We expect that 
there may appear phases with vortex lattice structures different from ordinary triangular 
Abrikosov lattice. 

The non-Abelian quantum turbulence phase where vortices are tangling and reconnecting 
here and there will be also possible if there is some mechanism which leads to an instability 
of lattice structures 



32J. It is interesting to investigate the property of this turbulence 
since we can expect behaviors which can be considerably different from atomic superfluids 
or helium superfluids. The energy spectra of turbulent fluids are determined by what kind 
of interaction transfers the energy from lower to higher momentum region. In a helium 
superfluid the Kelvin wave cascade and the phonon emission are considered as the dominant 
mechanism of energy transfer at the length scale which is smaller than the mean vortex 
distance. On the other hand, in the non-Abelian quantum turbulence, not only the Kelvin 
wave cascade and the emission of the U(1)b phonons from the vortex but also waves of the 
CP 2 orientational zero modes can contribute to the energy cascade. The propagation of CP 2 
waves can be affected by emission and absorption of virtual gluons. The dual Lagrangian 
obtained in this paper will provide a basis for the analysis of properties of the non-Abelian 
quantum turbulence. 

Finally the method we have used here can be applied to other types of non-Abelian vor- 
tices. In the models with gauged U(l) symmetry, non-Abelian vortices also exist which are 
called local non-Abelian vortices. In the literature, local non-Abelian vortices in relativistic 



field theories with (or without) supersymmetry have been studied extensively [23] (see [24 1 
for a review). In this case, vortices are coupled to only massive fields [33]. On the other 
hand, at sufficiently high density where U(1)a symmetry is effectively restored, global non- 



Abelian vortices 



26| become topologically stable. They are expected to interact with CFL 



mesons topologically. Via a dual transformation, CFL mesons are described by massless 
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non-Abelian two-form fields with antisymmetric tensor gauge invariance, which is known 
as the Freedman-Townsend model 1^. Unlike the case of semisuperfluid non-Abelian vor- 
tices, interaction between two non-Abelian global vortices was shown to depend on their 
<£P N ~ l orientations |27|]. This should be explained by the exchange of massless non-Abelian 
two-form fields. Extensions to these cases will be reported elsewhere. 
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Appendix A: Derivation of (1521) and (153ft 

Here we derive Eqs. (152 p and ( I53p . We make use of the property of Casimir operators. 
The generators of SU(iV) in the fundamental representation obey the following relation 

cnscnf = \^ k 3 -^^t (ai) 

First we show that (3> -1 )^t can be written as 

J = i 2 Tr (^ & ) - 2(1 - (3)^T a T b <p + 2(1 - - 7)0^00^} j (A2 ) 

with /3(r) = y^-p and j(r) = c+d^IJn) • us define the quantity 

X® b = ^T a T ai (p^T ai T a2 ^ ■ ■ ■ ^T an - x T h (\). (A3) 

By using the property of the Casimir operator ( 1A1I) and (fy<f> = 1 it can be shown that X n 
satisfy the following recurrence relation 

1 1/1 1 \"" 2 

Xf = -X? , ^T a (t)^T b <i). (A4) 

2 n_1 2N \2 2N J y ' 
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This equation can be solved to give 

n-l 



If = [^T a T b (j) - 2 {1 - (1 - I /NY' 1 } 0+^00^0] , (A5) 

with the initial condition X® b = <^T a T h ^>. The quantity * s defined by an expansion 

( *- i) --7^K^ + £(-^)" x " i )- (A6) 

Substituting Eq. (IA6j) into the above equation, now we perform the sum with respect to n 
and obtain the result (1A2|) . 

Then let us proceed to calculate {&~ x J) a = J 6 '^. The current for a = 

1, ■ ■ • , iV 2 — 1 can be written in terms of J 7 as Eq. ( 15TT) . Note that the following relations 
hold 

<p \T a T h ^T h <P = l - (l - 1) 0t T a 0; (A7) 



0tT a T fc 0Tr {T^(-l, 1)} = l0t T ^ ( _ 1; 1)0 _ _L t T>Tr ^ ( _ 1; 1} 

= iTr{T^(0, 1)}, 

0+rV.Tr {T b ^(-1, 1)} = 1)0 - ^ Tr ^(-1, 1) 

= 0, 



(A8) 



(A9) 



where we have used Eq. (lAip . 0^0 = 1 and the properties of J 7 . Applying these relations to 
the product of Eqs. ( 1A2I) and ( l5~Ti) leads to 

J); = -^70^0 {^0 + f|%o} + »*« Tr [T^(-l, /?)] , (A10) 

with a = ^r. By taking the complex conjugate, we have 

(J*~X = -^70^0 {<V? + ^|%o} - ^«Tr [T a ^(/3, -1)] . (All) 
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[28] U(1)b superfluid vortices found in [5|, |6j have winding number one while non-Abelian semisu- 
perfluid vortices have winding number 1/3, and consequently the former decay into three 
of the latter. Color magnetic flux tubes discussed in 6J, [7| have triple amount of fluxes of 
non-Abelian semisuperfluid vortices, and are topologically and dynamically unstable. 

[29] The time-dependent GL Lagrangian is valid when the temperature is close to the critical 
temperature and deviations from equilibrium are small. 

[30] The term Tr [d^ifi'ip — ^ d^tpj automatically vanishes since ip can be decomposed as ip = 
(A + p)1n + {x a + i( a )T a and the modes ( a are absorbed by gluons. 

[31] It is not possible to excite an on-shell gluon in terms of the CP N ~~ 1 effective theory. In order 
to describe the radiation of on-shell gluons from orientational modes, we have to include 
terms with higer derivatives in the CP^ -1 effective theory. However, we can still discuss the 
scattering of on-shell gluons by vortices in terms of the Lagrangian obtained in the paper, 
since there is no need to excite on-shell gluons in this case. 

[32] It has been shown in 



221 ] that, in the case of local non-Abelian vortices in the model with 
gauged U(1)b symmetry, collisions of two vortices always result in their reconnection even 
when they have different internal orientations. The same should hold for non-Abelian semisu- 
perfluid vortices because whether the reconnection occurs or not depends on only topology. 
[33] When the number of flavors is greater than the number of color, vortices are called semilocal. In 
this case, semilocal non-Abelian vortices are also coupled to massless Nambu-Goldstone bosons 
corresponding to spontaneously broken flavor symmetry [25[ as well as massive particles. 
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